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ABSTRACT 

In this paper, we review the observational properties 
of two accreting X-ray pulsars, the persistent source 
4U 1907+09 and the transient V 0332+53. Accreting X- 
ray pulsars are among the brightest sources in the X-ray 
sky and are frequently observed by Integral and other X- 
ray missions. Nevertheless they are still very enigmatic 
sources as fundamental questions on the X-ray produc- 
tion mechanism still remain largely unanswered. These 
questions are addressed by performing detailed temporal 
and spectral studies on several objects over a long time 
range. 

Of vital importance is the study of cyclotron lines as they 
provide the only direct link to the magnetic field of the 
pulsar. While some objects show cyclotron lines which 
are extremely stable with time and with pulse phase, in 
other objects the lines strongly depend on the pulse phase, 
the luminosity of the source, or both. 

Of special interest in any case are transient sources where 
it is possible to study the evolution of the cyclotron line 
and the source in general from the onset until the end of 
the outburst through a wide range of different luminosity 
states. 

We present an observational review of a transient and a 
persistent accreting X-ray pulsar observed with Integral. 

Key words: X-rays: stars - stars: magnetic fields - stars: 
flare - stars: pulsars: individual: 4U 1907+09 - stars: 
pulsars: individual: V0332+53 . 



1. INTRODUCTION 



Accreting X-ray pulsars are among the most prominent 
sources in the X-ray sky. Since accreting X-ray pul- 
sars are usually young objects, they are clustered along 
the galactic plane. Although already discovered in the 
first days of X-ray astronomy, many fundamental ques- 
tions remain unanswered to date. The X-ray production 



mechanism is only roughly understood, as a detailed un- 
derstanding would require a full magneto-hydrodynamic 
treatment, incorporating the presence of a magnetic field 
of ^10 12 G, and a full handling of general relativity. Es- 
pecially the behavior and shape of cyclotron resonant 
scattering features (CRSFs or cyclotron lines) are not 
very well understood and require detailed studying. But 
also only a very general picture of the overall accretion 
geometry and the accretion process itself exists. Analyz- 
ing the temporal behavior of accreting neutron stars along 
with the spectral properties allows to study these impor- 
tant questions. 

INTEGRAL is an excellent instrument to study these 
sources as INTEGRAL frequently observes the galactic 
center and the galactic plane and due to its large field of 
view. As a result, an enormous amount of data on accret- 
ing X-ray pulsars is obtained. In this paper we present 
some recent results on the analysis of two accreting X- 
ray pulsars: the detection of a complete torque rever- 
sal in 4U 1907+09 (3] and the evolution of the transient 
V 0332+53 over its recent outburst and especially the be- 
havior of the cyclotron line 11271 l33tl . 



2. 4U 1907+09 



The wind-accreting High Mass X-ray Binary system 
4U 1907+09 lfl3ll consists of a neutron star in an ec- 
centric (e = 0.28) 8.3753 d orbit around its companion. 
Based on a lower limit for the distance of 5 kpc, the X-ray 
luminosity is approximately 2 10 36 erg s _1 1151 . The stel- 
lar companion has been classified as a 08-09 la super- 
giant with an effective temperature of 30500 K, a radius 
of 26 Rq, a luminosity of 5 10 5 L Q , and a mass loss rate 
of 7 10 -6 M Q yr~* @. 

Similar to other accreting neutron stars, the X-ray con- 
tinuum of 4U 1907+09 can be described by a power-law 
spectrum with an exponential turnover at ~ 1 3 ke V 11251 . 
The spectrum is modified by strong photoelectric absorp- 
tion with a column Ah = 1.5-5.7 1 22 cm -2 ||7|] which 
is due to the dense stellar wind and is therefore strongly 
variable over the orbit. The column density is maximal 
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Figure 1. Light curve of4U 1907 +09, folded with the or- 
bital period of 8.3753 days [from 9]. Upper Panel: Rossi 
X-ray Timing Explorer All Sky Monitor (ASM) 2-12 keV. 
Lower Panel: IBIS 20-40 keV light curve. 



between the end of the primary X-ray flare and the start of 
the secondary flare (Fig. [T} of the orbital light curve JUJ]. 
At higher energies, the spectrum exhibits cyclotron reso- 
nant scattering features (CRSF) at ^19keV and ^40keV 
[USB]. 

With a pulse period of ~ 440 s, 4U 1907+09 is a slowly 
rotating neutron star j23ll . The neutron star has ex- 
hibited a steady linear spin down with an average of 
P pulse = +0.225 syr" 1 from P pulse = 437.5s in 1983 
to 440.76s in 1998 [16]. Recently, a decrease in P pu i se 
to ^0.115 syr -1 , about half the long term value was re- 
ported n. 



2.1. Timing Analysis 

The X-ray light curve of 4U 1907+09 clearly shows 
pulsations with a period of ^441 s. Additionally, 
4U 1907+09 also exhibits flares with a typical length of 
several hours 11231 [lql . Fritz et al. D9Q observed four such 
flares with count rate increases of at least 5u over its nor- 
mal level (Fig.[2]l. Three of them, in revolutions 187, 193, 
and 305, are associated with the main peak in the orbital 
light curve (Fig. [TJ while the flare observed in revolution 
188 is linked to the secondary peak. 

4U 1907+09 showed a steady spin down rate of P pu ise = 
+0.225 syr" 1 d. Baykal et al. 0] using RXTE-data 
determined a spin down rate 0.5 times lower than the 
long term value. Fritz et al. using INTEGRAL data, 
extend the tracking of the pulse period evolution to the 
years 2003-2005. Therefore these authors first determine 
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Figure 2. IBIS (ISGRI) 20^10 keV light curve of revolu- 
tion 193 [from 9]. The bin time is 200 s. The inset shows 
a close up of the flare, the bin time in this case is 40 s. 
X-axis numbers are hours. 
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Figure 3. Evolution of the 4U 1907+09 pulse period over 
the last 20 years [as shown by®. Diamonds indicate the 
results obtained in this work, other symbols indicate data 
from ft23\ plus-signs ], if^j asterixes ], ^25\ triangles ], i f/61 
squares], and £J crosses]. Note the clear deviation 
from the long term trend leading to a complete spin re- 
versal. The inset shows a closeup of the periods found by 
INTEGRAL Ml. 



the pulse period by first using epoch folding l22ll to get 
a starting period and then performing pulse phasing to 
get a higher accuracy. Fig. [3] shows the long-time his- 
tory of the period evolution based on all available data. 
Note that the recent periods obtained fromRXTE-data [2] 
show a clear deviation from the historic spin down trend 
of Pp U i S e = +0.225 syr -1 from IU5I1 . The periods ob- 
tained by Fritz et al. [9] from INTEGRAL -data (shown 
as diamonds) not only confirm this change of the trend, 
but show a complete trend reversal from the historic long 
term spin down trend to a spin up trend from MJD 53131 
onwards (see Fig. 0). 

Using the derived periods, Fritz et al. [9] obtained pulse 
profiles for all revolutions (Fig. |4|i by folding the light 
curves with the best respective periods. No change in the 
shape of the pulse profile is obvious, especially no change 
in the shape before, during, and after the reversal of P can 
be seen. 



2.2. Discussion 

In the INTEGRAL data, a spin up phase of 4U 1907+09 
is clearly detected after 20 years of spin down (Fig. [3]) at 
a rate of P pu i se = +0.225 syr -1 , with the available data 
being consistent with P = 0s _1 . In 2002 and 2003 first 
indications for a decrease in the magnitude of spin down 
were detected yj]. As shown in Sect. 12. II the new INTE- 
GRAL data show a torque reversal from ^MJD 53131 
onwards, with the source now exhibiting a spin up with a 
rate of P pu i se = -0.158 syr" 1 . So far, the INTEGRAL 
results appear to be consistent with P = s _1 (the la up- 
per limit for P during the spin up is —4 10~ 5 s _1 yr -2 ). 

The conventional interpretation of spin up in accreting 
X-ray pulsars with a strongly magnetized, disk-accreting 
neutron star is that the accretion disk is truncated by the 
magnetic field of the neutron star [see, e.g., [HI UH EH 
and references therein]. The transfer of angular momen- 
tum from the accreted matter onto the neutron star will 
lead to a torque onto the neutron star resulting in spin 
up. The general expectation of these models is that the 
Alfven radius r m , where the accretion flow couples to the 
magnetic field, is assumed to be close to the co-rotation 
radius, r co . 

This simple model, however, is unable to explain the long 
episodes of constant P. For 4U 1907+09 the magneto- 
spheric radius inferred from the cyclotron line measure- 
ments is r m ~ 2400km, while r co ~ 12000km jlUl . 
Moving r m out to r co would require a magnetic field of 
~10 14 G as for a magnetar, which is two orders of mag- 
nitude larger than the magnetic field deduced from the 
observed cyclotron lines. 

Fritz et al. @], however, then discuss that the a model pre- 
sented recently by Perna et al. B28I1 is able to explain these 
observations. This model is based on the fact that X-ray 
pulsars are oblique rotators, i.e., the neutron star's mag- 
netic field is tilted by an angle \ with respect to the axis 
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Figure 4. Pulse profiles for different revolutions in the 
20-40 keV band [from 9]. Panels a-e show pulse profiles 
obtained during the period of spin down, and panels f-o 
show spin up pulse profiles. 



of rotation of the neutron star. Assuming the accretion 
disk is situated in the neutron star's equatorial plane, for 
a ring of matter in the accretion disk, the magnetic field 
strength then depends on the azimuthal angle and thus 
the boundary of the magnetosphere is asymmetric. Such 
a configuration can lead to regions in the disk where the 
propeller effect is locally at work, while accretion from 
other regions is possible. This then results in a nonlin- 
ear dependence of the accretion luminosity from the M 
through the outer parts of the disk. Perna et al. 12811 point 
out that a nontrivial consequence is that for values of \ 
greater than a critical value, Xcnt, limit cycles are present, 
where cyclic torque reversal episodes are possible with- 
out a change in M. For typical parameters. Xcrit is be- 
tween ~25° and —45°. Perna et al. HH and Fritz et 
al. stress that the existence of the limit cycles de- 
pends only on \ an d on me pulsar's polar magnetic field, 
J5, and no variation of external parameters is required 
to trigger torque reversals, while the torque reversal of 
4U 1907+09 with no associated drop in luminosity and 
no change in the shape of the pulse profile is difficult to 
explain in conventional models. 



3. V 0332+53 



The recurrent X-ray transient V 0332+53 experienced an 
X-ray outburst from December 2004 to February 2005 
lfl9l l29ll . The outburst was predicted almost one year be- 
fore the actual outburst 111411 due to the optical brightening 
of the optical companion BQ Cam, a 08-09e star. 

The observation of this outburst of V 0332+53 b y R XTE 
revealed the presence of three cyclotron lines 141 1191 12911 . 
The fundamental CRSF has an energy of 28 keV, and the 
first and second harmonics are at 50keV and 71 keV, re- 
spectively. This corresponds to a magnetic field strength 
of 2.7 x 10 12 G Ella]- The spectrum is otherwise de- 
scribed by a power law modified by a high-energy cut-off 
at higher energies. 



3.1. Fluxes 

The light curve of V 0332+53 in ISGRI and JEM-X dur- 
ing the outburst is shown in Fig. [5] Mowlavi et al. i27tl 
describe the decline phase by an exponential decay up to 
MJD~53412, followed by a linear decrease. The expo- 
nential decay is less rapid at the ISGRI energies than at 
the JEM-X energies. From a simple estimation, the expo- 
nential decay time is about 30 days above 20keV (solid 
curves in Fig.|5}i) and about 20 days below 15 keV (solid 
curves in Fig. |5J)). The exponential decay is also con- 
firmed by the 2-12 keV light curve as observed by the All 
Sky Monitor (ASM) on board of RXTE (see Fig.^p). 

Reig et al. [30] discovered a very interesting behavior 
of the color-color diagram over the outburst (see Fig. [6]). 
They derived background-subtracted light curves corre- 
sponding to the energy ranges c\ = 4 — 7.5 keV, C2 — 




53350 53360 53370 53380 53390 53400 53410 53420 



t (MJD) 

Figure 5. a Light curve of V0332+53 over the decay of 
the outburst in different energy bands in Crab units [as 
shown bv \2%l . The solid curves are fits to the exponential 
(from MJD=53376 to 53412 with t=30 days) and linear 
(from MJD=53412 to 53422 ) b Same as a, but for JEM-X 
fluxes and with an exponential decay of t=20 days. The 
2-12 keV RXTE/ASM light curve is also shown, using a 
normalization factor of 77 cps for 1 Crab. 
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Figure 6. Color-Color Diagram of the outburst of 
V 0332+53 [from\3^ using RXTE-data. The X-ray col- 
ors are defined as soft color SC=7.5-10 keV/4-7.5 keV 
and hard color HC=15-30 keV/10-15 keV. 
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Figure 7. Evolution of the parameters of the fundamen- 
tal cyclotron line [fromWftl. a line energy, b width, and c 
depth. The vertical bars represent 90% uncertainty, while 
the horizontal bars indicate the duration of each obser- 
vation. 



7.5 - 10 keV, c 3 = 10 - 15 keV and c 4 = 15 - 30 keV 
and then defined the soft color (SC) as the ratio C2 /ci and 
the hard color (HC) as the ratio between C4/C3. 

These authors observed that the source was in a soft state 
during the peak of the outburst and as the count rate de- 
creased the source moved up and right in the color color 
diagram, i.e. it became harder. At the end of the outburst 
(regions E and F) the source moved again to the left, i.e., 
lower values of SC, while HC remained mostly constant. 
Therefore two spectral states or branches can be distin- 
guished: a soft state when the source is bright and a hard 
state when the source gets dimmer. 

The pattern in the color color diagram of V 0332+53 re- 
sembles that of low mass X-ray binary Z source [see e.g 
l34ll . In this picture, the soft branch would correspond to 
the normal branch and the hard branch to the horizontal 
branch in a Z source, however, as Reig et al. [30] point 
out, the flaring branch observed in Z source is missing 
in V 0332+53. During the decline of the outburst, the 
source moves without jumps through the Z shape in the 
color color diagram, taking about 80 days from the softest 
to hardest point [30]. 



3.2. Spectrum 



The spectra of accreting X-ray pulsars are a superposi- 
tion of contributions from the hotspots at both magnetic 
poles and the lower parts of the accretion columns. Since 
no self consistent model is available, various empirical 
models have to be used lfl8ll 
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Figure 8. a Spectrum and folded model for JEM-X 
(left) and ISGRI (right) using data of revolution 284 
[from \27\l : b residuals of the model with T=— 0.41 and 
E{ id=7.4 keV, and without any cyclotron line modeled; c 
one Gaussian is included at 29.3 keV; d a second Gaus- 
sian is added at 51 keV to model the first harmonic; e 
as before, but with the energy of the first line fixed at 
27.5 keV, the value found in revolution 273; f same as d, 
but with r fixed to the value of ^—0.17 found in revolution 
273. The strong residuals in the two bottom panels show 
that the parameters of the fundamental cyclotron line do 
indeed change. 
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Figure 9. Dependence of the energy of the energy of the 
fundamental cyclotron line of "V 0332+53 ffrom \33\l . Tri- 
angles represent results from INTEGRAL, and squares 
from RXTE. 



Mowlavi et al. [2jt] use the simple "cutoffpl" model to 
fit the combined ISGRI and JEM-X data, to which two 
cyclotron lines with a Gaussian optical depth profile JH 
have been added. A multiplicative constant is applied to 
allow for the different normalizations of the two instru- 
ments. Furthermore a systematic error of 2% is applied 
to account for the uncertainties in the response matrices 
of JEM-X and ISGRI. No photoelectric absorption is re- 
quired at low energies, compatible with the fact that the 
JEM-X data are analyzed only above 4keV. The same 
model is used for all observations to keep consistency be- 
tween the parameters from one observation to the next. 

The cyclotron line parameters for the fundamental line 
j27tl are displayed in Fig. [7] as a function of time. The 
spectrum with the fits and the residuals for a model with 
no, one, two, and three cyclotron lines using the data of 
revolution 284 is shown in Fig. [8] The energy of the fun- 
damental cyclotron line increases by 7% from the first to 
the last INTEGRAL observations, and the depth of the 
line increases by 16% while at the same time the width 
decreases by 12%. These changes are significant; accord- 
ing to Mowlavi et al. j27ll . the attempt to fit the spectrum 
in revolution 284 with a line at 27.5 keV results in strong 
residuals around the line (Fig. [8^), but also variations of 
continuum can be ruled out to explain the changes in the 
line parameters (Fig.JHf). 

Tsygankov et al. ll33ll then went one step further and in- 
stead plotting the energy of the fundamental cyclotron 
line versus time, these authors show the dependence of 
the cyclotron energy on the luminosity of the source, thus 
combining Figs. [7] and [5] into Fig. [9] A linear fit to this 



dependence results in the following relation: 

£ C yc,i = -O.IOL37 + 28.97keV, 

where L-yj is the luminosity in units of 10 37 erg s _1 . 
Assuming that for low luminosities the emission comes 
practically from the neutron star surface, Tsygankov et 
al. IB3I1 estimate the magnetic field on the surface of the 
neutron star for canonical values for the radius and the 
mass of the neutron star to B = 3.0 10 12 G. 



3.3. Discussion 

V 0332+53 is a highly magnetized pulsar with a magnetic 
field of about 3.0 x 10 12 G EH10, as derived from 
the energy of the fundamental cyclotron line. The X-ray 
outbursts from this system are observed with a time recur- 
rence of about 10 years. They are believed to be triggered 
by a massive mass ejection from the optical companion 
star, resulting in the optical brightening of the system 
in 2004 O. The ejected mass feeds an accretion disk 
around the pulsar, from which it is funneled by the mag- 
netic field towards the polar caps and falls on the surface 
of the neutron star. 

The analysis of the 2005 outburst provides new insights 
towards a better understanding of V 0332+53 and simi- 
lar transient sources within this basic picture. The X-ray 
light curve (Sect. 13. jj displays a characteristic exponen- 
tial decay with a time scale of 20 to 30 days, followed by 
a linear decrease. This remarkable change occurs approx- 
imately at the same moment when the trend in the color 
color diagram changes. This may provide some informa- 
tion on the accretion and especially the geometry. 

As discussed by Mowlavi et al. M27I1 . an exponential de- 
cay of the flux is observed also in dwarf novae and in soft 
X-ray transients on time scales from a few days for the 
former to several weeks for the latter j2ll l32ll . King & 
Ritter [17] show how illuminated disks in low-mass X- 
ray binary stars can produce such an exponential decay 
on the time scales observed for the soft X-ray transients. 
Subsequently, numerical calculations by Dubus et al. [8] 
showed, that indeed irradiated accretion disks produce an 
exponential decay of the X-ray luminosity. 

As Mowlavi et al. l27ll point out, it is tempting to make 
the parallel with the results observed in V 0332+53. 
The exponentially decreasing flux followed by a linear 
decrease is very similar to what is observed in disk- 
illuminated systems. Such a behavior as observed by 
Mowlavi et al. i27ll has never before been observed in 
a high mass accreting X-ray binary pulsar system. The 
X-ray flux is proportional to the mass accretion rate, and 
therefore the accretion rate would then also be propor- 
tional to the mass of the disk. An unknown change in 
the disk would then trigger the switch to the linear de- 
cay phase. Furthermore, Mowlavi et al. f27ll suggest 
that the two different decay times observed in V 0332+53 
are due to the presence of at least two regions contribut- 
ing to the X-ray continuum: one time scale governs the 



rate of mass flow onto the neutron star through the disk, 
leading to the exponential decay followed by a linear de- 
crease (see above); another time scale governs the spec- 
tral changes of the emission that may come from different 
heights within the accretion column ll27ll . 

The X-ray spectrum (Sect. 13.2b is well fitted by a cutoff 
power law modified by two cyclotron lines. The energy 
of the fundamental cyclotron line is shown to increase by 
20% from the peak of the outburst as measured by RXTE 
until the end of the outburst. At the same time it became 
narrower and deeper. This information would provide 
some insight into the resonance scattering region near the 
polar caps. The evolution of cyclotron line parameters 
with time has already been reported in the literature for 
several binary systems fl25ll . These authors attribute the 
change of the line energy to a change of the height of 
the scattering region above the neutron star, resulting in 
a height change ~300m for V 0332+53. Rreykenbohm 
et al. [20] report a variation of the cyclotron energy of 
GX 301-2 with the pulse phase of the pulsar, and attribute 
the change to different viewing angles of the accretion 
column where the line originates. 

For V 0332+53 Mowlavi et al. [27] interpret the in- 
crease of the energy as a variation of the scattering region 
characteristics. They assume that the region where the 
cyclotron resonance scattering takes place is extended. 
The region would then cover a range of magnetic field 
strengths resulting in a superposition of several narrow 
cyclotron lines, each originating at a specific height and 
therefore in a specific magnetic field strength, thus giv- 
ing rise to a broad line. As the accretion rate and hence 
the flux decreases with time, the cyclotron resonance re- 
gion shrinks and moves closer to the neutron star. At the 
end of the outburst only the region closest to the neutron 
star surface, where the magnetic field is strongest, would 
contribute to the formation of cyclotron lines. The re- 
sulting lines can therefore be expected to be observed at 
higher energies and being much narrower as is indeed ob- 
served in V 0332+53. Mowlavi et al. B27T1 estimate the 
total movement of the cyclotron line formation region to 
be of the order of 500 m, similar to 4U 01 15+63 where a 
change of ~1 km was observed 12611 . 
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